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Abstract 

Strutinsky's averaging(SA) method is applied to multiquark droplets 
to systematically extract the smooth part of the exact quantal energy 
and thereby the shell correction energies. It is shown within the bag 
model that the semi-phenomenological density of states expression 
given upto curvature order is almost equivalent to the SA method. 
A comparative study of the bag model and the relativistic harmonic 
oscillator potential for quarks is done to investigate the quark mass 
dependence of the finite-size effects. It is found that there is an im- 
portant difference between these two cases, which may be related to 
the presence/non-presence of the net spin-orbit effect. 

PACS numbers: 12.3.Mh, 12.39.Ba, 24.85.+p, 25.75. +r 







I. Introduction 



The most interesting prediction in recent times is the Witten's conjecture[0] 
that the strange quark matter(SQM) might be the absolute ground-state 
of hadronic matter. Such a possibihty has very significant astrophysical 
consequences 0. In addition, probable formation and detection of small 
lumps of SQM in relativistic heayy-ion collisions is being promoted as an 
unambiguous signature of quark-gluon plasma^]. In this context, stability 
of strangelets with baryon number A < 100 is of immense interest. In the de- 
termining the stability criteria of these small strangelets the finite-size effects 
such as surface and curvature play a significant role. 

Farhi and Jaffe]^ were the first to consider the surface effects. Later 
on, Mardor and Svetitsky|^ suggested that curvature effects are also im- 
portant. Recently, Madsen[^] has given a semi-phenomenological density of 
states(DOS) expression including terms upto curvature order for a system of 
quarks with mass rriq. He has shown within the bag model that this DOS 
expression reproduces quite well the exact shell model calculations irrespec- 
tive of the quark mass rriq. (We hereafter refer to this model as quark liquid 
drop(QLD) model.) Further, an intriguing aspect is that the finite-size con- 
tributions in the QLD model do not seem to show a converging trend. More 
specifically, in the case of massless quarks the surface energy is zero, whereas 
the curvature contribution is significant. And in the case of massive quarks, 
the surface energy contribution is non-zero while the curvature energy con- 
tribution can be zero or positive or negative depending on the value of quark 
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mass rriq. ( These features are clearly illustrated in Fig.l of Ref. ). Con- 
sequently, the question of the importance of still higher-order terms such as 
Gauss curvature naturally arises. In other words, is the agreement between 
the QLD and shell model calculations accidental? 

This question can be answered rather quantitatively with the help of 
the famous Strutinsky's averaging(SA) method|0. In the context of nu- 
clear physics, the SA method and the Wigner-Kirkwood(WK) expansion[§] 
have been shown analytically to be equivalent. Therefore, by comparing 
the QLD model and the SA method we can determine the importance/non- 
importance of higher-order terms in the semi-phenomenological DOS expres- 
sion of Madsen . In addition, the QLD model suggests that the finite-size 
effects depend on the value of the quark mass. As this can have very sig- 
nificant consequences in the study of quark-hadron phase transition 0, it 
would be interesting to know whether this finding is of general nature or 
is model-dependent. In this work, we mainly focus upon these two aspects 
of multiquark droplets. For this purpose, we consider the two potentials; 
the infinite square well(ISqW), namely a spherical cavity and the relativistic 
harmonic oscillator (RHO). The beauty of these two potentials is that both 
mimic quite adequately the asymptotic freedom and confinement properties 
of QCD. Moreover, they are analytically solvable. 

In section II, we briefly discuss the nature of the shell structure in both 
these potentials. The Strutinsky's averaging method appropriate for quark 
systems is presented in section III and the results obtained are discussed in 
section IV. Finally, we summarize our findings in section V. 
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II. Confining potential, quark mass and shell structure 



In regard to the shell structure in the ISqW and RHO potentials, at the 
outset the following observations can be made. 

• In the case of the spherical cavity, the eigenvalues obtained exhibit 
splitting of levels such as (^3/2,^1/2), (c?5/2, '^3/2) ■ ■ ■ due to the spin-orbit 
coupling. On the otherhand, in the harmonic oscillator the eigenvalues 
are j— independent and hence no net spin-orbit effect on the energy 
spectrum. Thus, as expected, the shell structure depends on the nature 
of the confining potential. 

• The other interesting observation is that the energy gap Aj between the 
/ = j + 1/2 and / = j — 1/2 levels decreases as the value of quark mass 
niq increases. This is illustrated in Fig.l where we have plotted the 



eigenvalues obtained by the bag boundary condition |jTO| taking iriq = 
0, 150, 450 MeV. The zero of the scale is taken at the Osi/2 level of 
massless case and therefore, what is plotted is (co'(k) — 2.0428) in units 
of hc/R. ( While solving the boundary condition for massive quarks 
we have taken R = 1 fm). The decrease in Aj can be clearly noted 
in the Op and Od shells. In addition, due to the dependence of Aj on 
rriq, rearrangement of levels occur at higher j values as compared to 
the massless quark spectrum. The degree of rearrangement ofcourse 



depends upon the magnitude of m^. In a recent work||TT[, it was found 
within the relativistic mean field theory of nuclei that the spin-orbit 
term is inversely proportional to the density dependent nucleon effective 



mass. In analogy, we may expect to decrease as the quark mass 
increases. Results displayed in Fig.l seem to illustrate this effect. On 
the otherhand, in the case of RHO potential, the shell structure remains 
more or less the same expect for a relative shift in energy as rUg is varied. 

Thus, already one can see that the nature of the confining potential and 
the quark mass have definite roles to play in determining the shell struc- 
ture and thereby the stability of small multiquark droplets. In the following 
section, we shall outline the SA method and apply it to a system of quarks. 

III. Strutinsky's averaging method for quark systems 

The basic concept of the SA method]^ is that the total energy E of a given 
quantal system can be separated into a smooth part E and an oscillating part 
AE, i.e. E = E + AE. For a multiquark system consisting of more than 
one flavor, AE pertaining to each flavor needs to be calculated separately. 
For the sake of simplicity, we consider here a system of one kind of quarks 
with mass nig confined either in a spherical cavity or in a harmonic oscillator 
potential. Then the exact single particle DOS is given by 



where g is the degeneracy factor and are the eigenvalues of the potential 
under consideration. The total energy of the quark system is then given as. 




(1) 



m 



m = occupied levels. 



(2) 



m 
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Now, to evaluate the smooth part of the total energy E, Strutinsky pro- 
posed a numerical averaging of the single particle spectrum g{e) by Gaussian 
smoothing funtions centered around and taken over a certain energy range 
7. Thus, the smoothed level density is given as 

M p 

m = ^ E expi-Ui)j:CMU^l (3) 

where Um = (e ~" ^m)/l, M is the number of levels taken into consideration 
and p = 6|g. The coefficients Q = {-if/^ /[T{i/2)\] for even i values, and 
Ci = for odd i values. Hi{Um) are the Hermite polynomials. The average 
number of particles N and E are then defined as, 

N = f'g{e)de, 

E = f\g{e)de. (4) 

Having obtained the exact value E and the smooth part E, we can deter- 
mine the shell correction energies using the definition AE = E — E. In the 
following, we apply this method to the RHO and ISqW potentials. 

A. Relativistic harmonic oscillator potential 

Considering the quarks with current quark mass nig to be confined by a 
scalar plus vector harmonic confinement of the type |(1 + '-fo){mq + C2r^), 
the corresponding Dirac equation can be solved analytially[|l^]. And the 
eigenvalues are determined by the equation 

(e^-m,)(e^-mJ)=4C2(A + ^)^ (5) 
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with = 2n+l. In the non-relativistic hmit, eN = (^+3/2)y {2C2/mq)+mq 

2 /3 

and in the hmit of small mass, e^v — [2y/C2{N + 3/2)] + rrig/S. In analogy 
with the non-relativistic case, we take the potential strength parameter C2 = 
C2o^~^^^, where A = Nq/3 and Nq is the number of quarks. The quantity 
C20 can be related to the ground-state energy of a 3 — g system as C20 = 
{l/9){Mb/3). In the case of massless quarks, = 1085.5 MeV and C20 = 
(173.95)3 MeV^ and for massive quarks, Mb = 1672.45 MeV and C20 = 
(268.01)3 MeV^. 

The next important step is to choose an appropriate value of 7, and check 
that the smoothed part does not depend upon 7 value. This can be done 
by using the so-called plateau condition |jl^. In our study, we choose the 7— 



parameter to be 7 = 706*20^ A"^/^^ Then, the value of AE corresponding to 
the plateau region pertaining to each value of rriq and A is taken to be the 
'physical' shell correction energy which then determines the value of E. Fur- 
ther, it was found that the prominence of the plateau and its length depends 
upon the number of quarks Nq. As A decreases, it becomes somewhat tricky 
to fix the value of AE. The indeterminacy is however small. 

It would be interesting to compare the so-determined values of E with 
those obtained using a Wigner-Kirkwood expansion and study the impor- 
tance of higher-order finite-size effects such as Gauss curvature. The expres- 
sions for the quark number and the total energy using the WK expansion can 
be obtained analytically|[l5|]. We state below the final expressions so-obtained 
upto 0{h 
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rx^ + — x + 1- 1 + x 



-1/2 



-1/2 



,(6) 



where x = {rj — mq)/mq and is fixed by the number equation. Before, we 
present our results in this case, we would like to discuss the SA method in 
the case of ISqW potential. 



B. Infinite square well potential 



Considering Nq number of quarks with mass rUq to be confined in a spher- 
ical cavity of size R, the boundary condition |T0[ to be solved for determining 
the eigenvalues is 

ui + fl 

where a;^ = tu^ + jj? with /i = rriqR, and 1 = 1 — sgn^K). Then, the total 
energy of the system using the bag model picture is, 

E=^Y.-^. + \'kR'B, (8) 
R . 3 

where B is the bag energy density. The equilibrium radius of the system is 
then determined by the saturation condition, dE/dR \r^= 0. 

In this case, the single particle DOS g{uj) = — ^i) is smoothed as in 
Eq.(3) and is given as, 

M p 

g{uj) = exp{-Ul)J2aH,{Um), (9) 

7wV ^ m=l i=l 



with Um = {u! — u!m)/ju>- Then the smoothed energy E can be obtained using 
the equation, 

E = ^K{R) + I-kR^b, (10) 

where E^ = ijjg{uj)duj and A^^ is determined by the number equation. 
Here, the 7^, is chosen as; 7^, = ^qA~^/^. 

With this parametrisation, we have checked the plateau condition in this 
case also for each value of ruq and A while calculating E. 

IV. Results and discussions 

Having demonstrated a reliable way of extracting the smooth part from 
the exact quantal energy, we shall presently compare this SA method with 
the WK one in the case of RHO potential, and the QLD model in the case 
of ISqW potential. 

In Table I, we have presented the energies obtained in the case of RHO 
potential for four values of baryon number A using the SA method [Eqs.(3-5)] 
and the WK method[Eq.(6)] and then are compared with exact ones[Eq.(2)]. 
The difference between E and E^k is about 10 MeV for massless quarks and 
about 17 MeV for the massive ones. This may be attributed to the presence 
of higher-order [0(/i^)] correction terms. Notwithstanding this, we might say 
that the agreement is quite good. This then establishes the goodness of the 
WK expansion for further study of the quark mass dependence of surface 
and curvature energies. 

In Table II, the results obtained in the case of the ISqW potential is given. 
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We have taken = 1 MeV and 5^/^=145 MeV. One can immediately see 
that Eqld agrees quite well with that of the SA method. Similar degree of 
agreement is also obtained for more massive quarks. This then demonstrates 
that the DOS expression upto curvature order of Madsen is quite adequate. 

Consequently, an important question arises; Does the finite-size effects 
greatly depend on the mass of the quark as suggested by the QLD model? 
To answer this, we need to understand the following aspects of the liquid 
drop model(LDM) expansion of energy [ll6|, [T7| . The surface and curvature 
energy coefficients contributes respectively to orders of A^^^ and A^^^ in the 
LDM expansion of total energy: 

E = ayA + ttsA'^/^ + OcvA^/^ + ■ ■ • . (11) 

To these orders, there is a part coming purely from the surface and curva- 
ture terms of the DOS expression. But, there is also some extra contributions 
arising out of the finite size effect on the Fermi momentum kp- ( This can be 
noted from Eq.(4) of Ref. ). Because of this, the volume term in the DOS 
expression contributes towards 0{A^^^), 0{A^^^) ■ ■ -, and similarly the sur- 
face term contributes towards 0(^4^/^) and 0(^4'^), and so on. Therefore, the 
vanishing of curvature term in the DOS expression does not necessarily mean 
that curvature energy coefficient is zero, but ofcourse its value is expected 
to be small. In view of this, we would like to extract the effective surface 
and curvature coefficients by making a least-squares fit to the total energies 
Eqld calulated using three values of m^. To do so, we used the QLD model 
with B^^^ = 145 MeV and have obtained energies for 101 number of droplets 
with baryon number A such that Q < A < 1000. The results so-obtained for 
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the energy coefficients are given in Table III for niq = 0, 150 & 450 MeV. 

It can be seen that the volume ay coefficients obtained here agree with the 
bulk limits shown in Fig.l of Ref.[^ for all the three values of m^. Further, 
in the case of massless quarks the surface energy coefficient is indeed zero, 
and is about 220 MeV. As the quark mass is increased to 150 MeV, 
the curvature coefficient has decreased as compared to the massless case 
illustrating the fact that the curvature term in the DOS expression is nearly 
zero for = 150 MeV. Similarly, the value of is now non-zero due to 
the presence of a surface term in the DOS expression. With further increase 
in trLg, the value of increases, while Ocv remains more or less the same. 
Hence, although the curvature contribution in the DOS expression becomes 
negative as rriq increases, the effective curvature energy remains positive. 
Further, except for the massless case there is a converging trend in the LDM 
expansion of energy. It must be said here that the estimates obtained for 
higher-order coefficients such as Gauss curvature {0{A^)) are not reliable as 
such terms were not included in the DOS expression. We have given them 
just to show that the three principal terms a^, and stabilises with 
respect to the number of parameters in the ffi. Thus, there is a systematic 
dependence of surface and curvature energies on the mass of the quark. It 
must be stressed that the values obtained here are dependent on the bag 
constant value. Now, we are curious to see how far these findings are true in 
the case of RHO potential. 

For this purpose, we repeated the same exercise taking the WK expansion 
of number and energy (Eq. (6)) for the same set of A and m^. Here also we 
have checked the stability of the three energy coefficients with respect to the 
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number of parameters in the fit. The results so obtained are given in Table 
IV. It can be seen that the surface coefficient is independent of the mass of 
the quark. And the curvature coefficient increases from 150 MeV to 220 MeV 
as ruq varies from to 150 MeV. With further increase in m^, there is a slight 
decrease in a^v In the case of massless quarks, it is somewhat disturbing to 
find the volume coefficient a„ less than 930 McV. This can be rectified, as 
done in the case of bag models in choosing B, by appropriately choosing the 
potential parameter C20 so that the ud matter is unbound against nuclear 
matter. But, we feel our findings regarding the quark mass dependence of 
the finite-size effects shall remain unaffected as the surface and curvature 
terms are dependent only on the difference between the total energy and the 
bulk value. 

Further, the surface energy coefficient remaining zero irrespective of the 
value of the quark mass needs more critical examination as it can have signif- 
icant consequences for the phase transition studies. It may be recalled here 
that in the case of ISqW potential, the shell structure depends on the quark 
mass through the spin-orbit effect; whereas, in the case of RHO potential 
there is no such effect. Is this the underlying reason for the weak dependence 
of Us and Ucy on rUq in the case of RHO potential? May be. 

Thus, we have shown that the shell structure as well as the quark mass 
dependence of the finite-size effects are dependent on the nature of the con- 
fining potential. 

VI. Summary 
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In summary, we have applied the Strutinsky's averaging(SA) method to 
multiquark droplets to systematically extract the smooth part of the exact 
quantal energy, and thereby the shell correction energies. It is shown in the 
case of bag model picture that the DOS expression given upto curvature or- 
der reproduces quite well the smoothed energies obtained by the SA method. 
Similarly, we found in the case of relativistic harmonic oscillator(RHO) po- 
tential, the Wigner-Kirkwood expansion with terms upto 0{h^) reproduces 
well the smoothed values of energy. 

Having established the goodness of the asymptotic expansions in both 
the cases, we then made a comparative study of the two potentials, namely 
the spherical cavity and the relativistic harmonic potential, in regard to the 
quark mass dependence of the finite-size effects. 

It was found in the case of the RHO potential that in contrast to the bag 
model picture, the surface and curvature coefficients are weakly dependent 
on the quark mass rriq. Further, the surface energy contribution is almost 
zero irrespective of the value of rUg. These differences may be traced back to 
the difference in the shell structure between these two potentials due to the 
presence/non-presence of net spin-orbit effect. 

Acknowledgements: Useful discussions with Professor J. Madsen and Dr. 
M.G. Mustafa are gratefully acknowledged. 
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TABLE CAPTIONS 



Table I : Relativistic harmonic oscillator potential. Total energies calu- 
lated using the Wigner-Kirkwood(WK) expansion and the Strutinsky's aver- 
aging(SA) method are compared with the exact quantal ones for four values 
of baryon number A. The first set corresponds to the massless quarks, and 
the next one corresponds to quark mass = 150 MeV. 

Table II : Infinite square well potential. Total energies calulated using 
the quark hquid drop(QLD) and the Strutinsky's averaging(SA) method are 
compared with the exact quantal ones for four values of baryon number A 
taking the quark mass mq=l MeV. 

Table III : Infinite square well potential. Values of the energy coefficients, 
volume a^, surface and curvature Ccv, obtained by making a least-squares 
fit to the total energies calculated for 101 values of the baryon number A in 
the range 6-1000 using the quark hquid model. The three sets corresponds 
to quark mass niq — 0, 150, & 450 MeV. 

Table IV : Relativistic harmonic oscillator potential. Values of the energy 
coefficients, volume a^, surface Us and curvature Ocv, obtained by making a 
least-squares fit to the total energies calculated for 101 values of the baryon 
number A in the range 6-1000 using the Wigner-Kirkwood expansion method. 
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The three sets corresponds to quark mass = 0, 150, & 450 MeV. 
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FIGURE CAPTIONS 



FIG.l : Single particle levels obtained by solving the bag boundary 
condition are shown relative to the Osi/2 level of the massless quark case in 
units of he/ R. Index = 1, 2, & 3 corresponds to the quark mass niq—O, 150 
and 450 MeV respectively. 
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Table I 



A 






K 


6 


5556.9 


5466.2 


5476.8 


20 


17597.3 


17713.5 


17723.6 


40 


34979.6 


35126.0 


35134.1 


70 


61013.4 


61189.4 


61196.4 


6 


9794.6 


9653.1 


9670.0 


20 


31275.6 


31449.0 


31460.7 


40 


62239.8 


62456.6 


62467.2 


70 


108629.3 


108883.9 


108899.7 



Table II 



A 


E 


EsA 


Eqld 


6 


6924.4 


6948.9 


6965.4 


20 


22580.2 


22429.0 


22437.7 


40 


44581.0 


44406.5 


44413.6 


70 


77231.6 


77291.1 


77296.4 
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Table III 



4t Pa ram 




r 








3 


1090.3 


1.258 


218.0 






4 


1090.4 


-0.5139 


225.3 


-8.45 




5 


1090.4 


-0.0296 


222.3 


-0.838 


-6.39 


3 


1212.4 


265.0 


91.8 






4 


1211.8 


273.7 


55.9 


41.4 




5 


1211.9 


272.8 


61.4 


27.7 


11.5 


3 


1837.8 


303.4 


78.0 






4 


1837.5 


307.5 


61.0 


19.6 




5 


1837.5 


306.8 


65.3 


8.78 


9.10 
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Table IV 



4t Pavam 










Uy 


3 


865.4 


-2.68 


161.1 






4 


865.1 


1.19 


145.0 


18.5 




5 


865.2 


0.015 


152.4 


-0.015 


15.5 


3 


1542.9 


-4.0 


233.8 






4 


1542.5 


1.79 


209.8 


27.7 




5 


1542.6 


-0.096 


221.7 


-1.93 


24.9 


3 


2175.7 


-3.56 


194.9 






4 


2175.4 


1.55 


173.8 


24.4 




5 


2175.4 


0.053 


183.2 


0.86 


19.8 
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